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Abstract. A computational method to obtain optical conductivities from first principles is presented. It exploits a relation 
between the conductivity and the complex dielectric function, which is constructed from the full electronic band structure 
within the random-phase approximation. In contrast to the Drude model, no empirical parameters are used. As interband 
transitions as well as local-field effects are properly included, the calculated spectra are valid over a wide frequency range. As 
an illustration I present quantitative results for selected simple metals, noble metals, and ferromagnetic transition metals. The 
implementation is based on the full-potential linearized augmented-plane-wave method. 
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INTRODUCTION 



THEORETICAL FRAMEWORK 



The optical conductivity a yields a wealth of information 
about the electrical conductance properties of a material 
and its reaction to applied electromagnetic fields in linear 
response. It is formally defined by the ratio 



j(q,co) = c7(q,co)E(q,co) 



(1) 



between the total electrical field E inside the material 
and the induced current density j. As all quantities are 
taken as macroscopic averages over many unit cells, Eq. 
(1) is only valid for small wave vectors q, otherwise the 
variations at the atomic scale must be explicitly included. 
In general, the conductivity is a tensor, but it reduces to 
a scalar for systems with cubic symmetry. 

The conductivity can be properly expressed within a 
quantum-mechanical context through the Kubo formula 
[1], but most actual studies of the optical and transport 
properties of solids are based on simpler schemes, such 
as the classical Drude model [2]. While its simplicity 
makes the Drude model very attractive, there are also 
clear disadvantages. In particular, it relies on parameters, 
especially the phenomenological relaxation time, that are 
difficult to establish unambiguously, because the optimal 
values may depend on the frequency range and on the 
specific optical function under consideration. Besides, 
the original Drude model ignores interband transitions 
in real inhomogeneous solids, which akeady dominate 
in the visible range, so that extensions with even more 
empirical parameters become necessary in this case. 

Here I describe a parameter-free ab initio method for 
calculating the optical conductivity of metals over a wide 
frequency range, which takes interband transitions and 
local-field effects in inhomogeneous solids into account. 
Quantitative results are shown for selected simple metals, 
noble metals and transition metals from the infrared to 
the ultraviolet region of the spectrum. 



Due to the small momentum transfer in optical absorp- 
tion it is usually sufficient to consider the limit q — >^ 0. 
The basic equation used in this work is the relation 



eu{(o) = 1 



Ani 



lima(q,«) 

0) q^O 



(2) 



between the conductivity in the long-wave-length limit 
and the macroscopic dielectric function, which follows 
from Maxwell's equations. The latter is defined as [3, 4] 

eu{(o) = lim — j-i (3) 

q^oe-i(q,a)) 

in terms of the matrix inverse of the microscopic dielec- 
tric function projected onto a plane-wave basis 

e-\(i,(o)^ j e-\ry-(o)e-"^-^'-^U\d\' . (4) 

However, actual applications often replace the left-hand 
side of (2) by limq^o£(q: which is computationally 
simpler, because it requires only a single element instead 
of the entire matrix, but neglects local-field effects. 

The microscopic dielectric function represents the link 
to the electronic structure of the material and can be 
rigorously derived in time-dependent density-functional 
theory [5, 6]. Here I use the random-phase approximation 

e(r,r';«) = 5(r-r')- / v(r-r")Zo(r",r';a))c/V" 



(5) 

with the Coulomb potential v(r — r') — l/|r — r'| and the 
susceptibility (density correlation function) of the non- 
interacting auxiliary Kohn-Sham electrons 

Xa{r,r; co) = Hm F L L (/"'k'<^ " f"^<^) 

' a n.kn'.k' 

^ 9«ka (r) 9„\'a (r) ^n'VJa (rQ 9n*ka i^') 
(0 - e„ka + £„'k'a + 
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The sums ran over all single-particle states with band 
number n, wave vector k, and spin G, which are obtained 
from a self-consistent solution of the Kohn-Sham equa- 
tion [7] of stationary density-functional theory [8]. The 
occupation numbers /„ka obey the Fermi distribution. 

While the method itself can be formulated in a very 
general way, it should be pointed out that the present im- 
plementation has some restrictions: The random-phase 
approximation (5) only includes the scattering between 
the electrons and the influence of the static crystal po- 
tential on the electronic stracture. This is appropriate for 
most optical experiments, but phonon scattering cannot 
be neglected in studies of the electrical resistivity. In such 
cases the effects of dynamical lattice vibrations on the 
dielectric function must be exphcitly taken into account 
[9]. Furthermore, the construction (6) assumes a perfect 
periodic crystal, where the Bloch wave vector k is a good 
quantum number, without defects or disorder. 

RESULTS 

For the calculations 1 employ the full-potential linearized 
augmented-plane-wave method. Exchange and correla- 
tion are treated in the local-density approximation [7]. 
The Kohn-Sham orbitals are expanded in plane waves 
with reciprocal lattice vectors up to Gmax = 4 Bohr" ^ 
in the interstitial region and up to angular momentum 
'max = 10 inside the muffin-tin spheres. I include second 
energy derivatives of the muffin-tin functions in order to 
improve the description of higher unoccupied states [10]. 
The mixed product basis for the response functions and 
the Coulomb matrix is constructed with Gf^ax = 5 Bohr^ ' 
and Lmax = 4. Integrations over the Brillouin zone are 
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FIGURE 2. Optical conductivity of copper compared to ex- 
perimental data from Ref. [13] (circles). 
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FIGURE 3. Optical conductivity of silver compared to ex- 
perimental data from Refs. [13] (circles) and [14] (squares). 
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FIGURE 1. Optical conductivity of aluminum compared to 
experimental data from Ref [12] (circles). In all figures the 
solid and dashed line indicate results with and without local- 
field effects, respectively. The inset shows the energy-loss func- 
tion, which is dominated by the plasmon resonance near 15 eV. 



performed with 32x32x32 mesh points, the singularity 
of the Coulomb matrix at the zone center q = is treated 
analytically [11]. AU spectra are calculated at the exper- 
imental lattice constants and at zero temperature. 

In the following 1 present results for selected metallic 
systems. As a first example. Fig. 1 displays the real part 
of the optical conductivity of a simple metal, aluminum, 
whose electronic structure resembles that of the homo- 
geneous electron gas. Nevertheless, interband transitions 
play an important role in the visible range and give rise 
to the resonances at 0.5 eV and 1 .5 eV. Both are predicted 
by this theoretical approach in very good agreement with 
experimental data from ellipsometry measurements [12]. 
Local-field effects have only a very small influence. The 
inset displays another independent optical quantity, the 
so-called energy-loss function Im[— l/eM(ft>)], which is 
dominated by the plasmon resonance near 15 eV. 

In comparison, the spectral features of noble metals 
exhibit more stracture, as shown in Fig. 2 for copper 
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FIGURE 4. Optical conductivity of iron compared to exper- 
imental data from Ref. [15] (circles). 
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FIGURE 5. Optical conductivity of nickel compared to ex- 
perimental data from Ref [15] (circles). 



and in Fig. 3 for silver. In these cases transitions from 
the occupied d bands dominate in the considered energy 
window and give rise to pronounced peaks in the visible 
to ultraviolet region when threshold values for transitions 
are reached. All the features measured by experimental 
spectroscopies [13, 14] are again very well reproduced. 

As shown in Fig. 4 for iron and in Fig. 5 for nickel, the 
spectra of transition metals are even more complex due to 
their comphcated electronic structure. In these materials 
interband transitions between the partially occupied and 
weakly dispersive d bands near the Fermi level give rise 
to a prominent peak structure at lower frequencies than 
in the noble metals. Furthermore, the exchange splitting 
in ferromagnets effectively doubles the band structure, 
leading to distinct optical transitions in the majority and 
minority spin channels that create additional features in 
the spectra. The theoretical curves are in good agreement 
with optical measurements [15], but they have sharper 
features, as no artificial broadening is appUed. 



CONCLUSIONS 

I have described a computational method for determining 
optical conductivities without any empirical parameters 
from first principles. Interband transitions and local-field 
effects are fully taken into account. Quantitative results 
for selected simple metals, noble metals, and transition 
metals confirm that practical calculations are feasible and 
that the method yields results in good quantitative agree- 
ment with experiments for materials with very different 
electronic characteristics. The essential spectral features 
are accurately described over a wide frequency range. 
Other optical quantities, such as the refractive index, the 
energy-loss function, or the reflectivity, can be derived in 
the same fashion from the complex dielectric function. 
Although the large numerical effort precludes a straight- 
forward apphcation to microscale device structures, such 
atomistic simulations of materials properties can be used 
to determine parameters for simplified models like the 
Drude model and its various extensions. 
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